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Abstract 
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in a good agreement with those found recently by exploiting the five-loop 
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1 Introduction 



The study of the critical properties of magnetic phase transitions in three 
dimensional (3D) cubic crystal is a problem attracting theoretical efforts over 
more than 25 years. By using the lower-order renormalization-group (RG) 
approach, Wilson and Fisher [p]], Aharony 0, and Ketley and Wallace ||. f| 
showed that in the critical region the fluctuation instability of continuous 
phase transitions may be observed, and that it may lead to the isotropization 
of the system with a cubic anisotropy. This fact gave rise to the question of 
what regime of the critical behaviour is actually realized in 3D cubic crystal 
with N = 3. It was soon understood that the calculation of the critical 
dimensionality N c of the order parameter is the crucial point in studying 
critical phenomena in a cubic crystal. Indeed, the critical value N c separates 
two different regimes of critical behaviour of the system. For N > N c the 
cubic rather than the isotropic fixed point is stable in 3D. At N = N c the 
points interchange their stability so that for N < N c the stable fixed point is 
the isotropic one. However, attempts to evaluate the critical dimensionality 
resulted in dramatically different estimates. 

In fact, the one-loop RG analysis of the stability matrix eigenvalues of 
the cubic and isotropic fixed points as well as some symmetry arguments 
(see section 2) predicts that N c should lie between 2 and 4. Many years 
ago the three-loop expansion for N c as a power series in e was obtained [[|. 
Summation of that short series at e = 1 (D = 3) by means of the Pade 
approximant [1/1] yielded the value N c = 3.128 |5J, while making use of 
the Pade-Borel resummation method results in the estimate N c = 3. In 
contrast to this, the value N c = 2.3 has been found on the basis of the 
variational modification of the Wilson recursion relation method in . Later, 
however, Newman and Riedel showed by decoupling the infinite system of 
the recursion relations for the scaling fields and then solving them that for 
D = 3 N c ~ 3.4 @. At the same time, the classical technique of high- 
temperature expansions, under some circumstances, allowed one to establish 
that for N = 3 the isotropic critical asymptotics in a cubic crystal is unstable 
[§], thus implying N c < 3. Ten years ago the analysis of the critical behaviour 
of the (mn)-component field model, which has a good number of interesting 
applications to the phase transitions in real substances, has been carried out 
within the three-loop RG approach in three-dimensions. The calculation of 
the stability matrix eigenvalues for the cubic model (m = 1, n = 3) provided 
the stability of the cubic fixed point in 3D, and the critical dimensionality 
turned out to be equal to 2.91 [f|. In agreement with this, the estimate 
N c = 2.9 was given in |I0|. More recently, Kleinert and Schulte-Frohlinde 



calculated the RG functions for the cubic model in (4 — e) dimensions up to 



five-loop order [11 1. Summation of the critical dimensionality expansion with 



the help of the Pade approximant [2/2] gave the estimate N c = 2.958 [[12 
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The cubic fixed point eigenvalues found by means of a simple resummation 
algorithm of the Borel type, accounting for the large-order behaviour of the 



(3-i unctions when the anisotropy parameter is very small [I3|, indicated that 



the cubic point is stable in 3D Finally, in the recent work Ref. [HJ by 
using finite size scaling techniques and the high precision Monte Carlo (MC) 
simulation it has been suggested that N c coincides with three exactly. Such 
strong scattering in the estimates of N c motivated us to study this problem 
with particular care. Calculation of the critical dimensionality as well as the 
eigenvalue exponents for the cubic and isotropic fixed points by exploiting the 
higher-order RG approach in three dimensions and generalized Pade-Borel- 
Leroy (PBL) resummation technique is the main goal of the paper. As will 
be shown, our estimates for N c and eigenvalues are in excellent agreement 
with recent results [15] obtained on the basis of the five-loop e-expansions. 



The layout of the paper is as follows. In the next section the model 
Hamiltonian is introduced and the massive field-theoretical RG procedure in 
fixed dimensions is formulated. The perturbative expansions for /3-functions 
for generic N are then deduced up to the four-loop order. In section 3 the 
structure of the RG flows of the model are investigated and the fixed point lo- 
cations are calculated using the generalized PBL resummation method. The 
eigenvalue exponents of the most intriguing 0(iV)-symmetric and cubic fixed 
points are evaluated for the physically significant case N = 3 and the stability 
problem is solved. The numerical estimate of the critical dimensionality N c , 
at which the topology of the flow diagram changes, is obtained by resumming 
both the four-loop RG expansions for the /3-functions in 3D and the five-loop 
e-expansion for N c at e — 1. In the conclusion the results of the investigation 
are discussed, along with the predictions and numerical estimates obtained 
earlier on the basis of the same or other theoretical approaches. 



2 The model, RG procedure and /3-functions 

We start from the fluctuation Hamiltonian 

( 1 1 

H = J d d x[-{m 2 <p 2 + d^id^i) + — (u G] jkl + v G 2 jkl )^ i ip j if k ipi] , (1) 

where tpi, % — 1, . . . , N, is the real vector order parameter field in fixed d and 
771q is the linear measure of the temperature, uq and vq denote the "bare" 
coupling constants. The symmetrized tensors associated with isotropic and 
cubic interactions are 

Gljki = ^(fiijfiki + 8ik8ji + 8u$kj) , G 2 ijkl = SijSikdu (2) 
respectively. 
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The model (|I]) has a number of interesting applications to the phase tran- 
sitions in three-dimensional simple and complicated systems. Indeed, when 
N = 1 Hamiltonian ([]]) describes the critical phenomena in pure spin system 
(Ising model), while for N = 2 it corresponds to the anisotropic XY model 
describing structural phase transitions in ferroelectrics as ordering the two- 
component alloys [Lj]]. The magnetic and structural phase transitions in a 
cubic crystal are governed by model (|IJ) as iV = 3. In the replica limit iV — > 
Hamiltonian ([J) is known to determine the critical properties of weakly disor- 
dered quenched systems undergoing second-order phase transitions [0 with 
a specific set of critical exponents |L9) . Finally, the case N — > oo corresponds 
to the Ising model with equilibrium magnetic impurities [20]. In this limit the 
Ising critical exponent of specific heat a changes its sign and takes the Fisher 
renormalization [21| together with v and 7: a —>■ —a/(l — a), v — > u/(l — a), 
7 -> 7/(1 -oc). 

To calculate the /3-functions normalizing conditions must be imposed on 
renormalized one-particle irreducible inverse Green functions 

} (P5 m,u,v; A; d) = Z^ /2 T iN) (p; m , tt , A; d) 

given by corresponding Feynman diagrams, A is the ultraviolet momentum 



cut-off. Within the massive field-theoretical RG scheme |E| at zero external 
momenta and at the limit A — > 00 they are normalized in a conventional way 



t r\p, ~V\ m,u,v;d) 
-^ t r\p, -p;m,u,v;d) 



ruR{{Pi}]m,u,v;d) 

rv R ({piY,rn,u,v-,d) 



p=0 
p=0 



{Pi}=0 



m 
1 



4— 

m u 



4—d 

m v 



(3) 



where m, u and t> are the renormalized mass and dimensionless coupling 
constants. The vertices are connected with the vertex function 

without external lines normalized at zero external momenta 



From equations (^]) the expansions for the renormalization constants Z^, Z u 
and Z v may be obtained 



Z" 1 



af?r (2) (p, -p; m ,M , ^o) 
^ 4) (0;mo,« ,t) ) , 



iiO 
«0 



(4) 



These constants relate the "bare" mass m and coupling constants Uq and t>o 
of the initial Hamiltonian (H) to the corresponding physical parameters 

ml + 5ml = m 2 , m = Z ip T^ 2 \p, -p;m ,u ,v ), 
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4— rf 



4-d Zv 

—u, v = m —v. 



(5) 



With relations 
the formulae 



taken into account, the /3-functions can be calculated via 

d In up o _i_ <91n up o 

dv Pv 
a In vp a _i_ d In vp g 



dv 



(d-4), 
(d-4), 



(6) 



where /3„ = 



9<? 



5- = 



f c*| lnm| ; 

For each Feynman graph contributing to the RG functions the corre- 



sponding contractions are computed by the algorithm developed in Ref. [24 



The combinatorial factors as well as the integral values are known from Ref. 



25(1 • After some work we obtain the four- loop expansions for the /^-functions 
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in three dimensions: 



0u 



U 



1 



u 



N + 8 



v + 



1 



(N + 



+ 44.444444 u v + 10.222222 v' 



3 (2.024691 N + 9.382716) u 2 
1 



(N 



3 (0.449648 N 2 



+ 18.313459 N + 66.546806) u 3 + 3 (6.646878 N + 164.613849) u 2 v 
+ 3 (0.621889 N + 100.955929) u v 2 + 65.937285 v 3 
1 



4 — — -(0.155646 N 3 - 35.820204 N 2 - 602.521231 N 

- 1832.206732) u 4 - 3 (1.352882 N 2 - 182.073890 N 

- 2064.170701) u 3 v + 3 (27.250336 N + 2110.408809) u 2 v 2 

+ 9 (1.291017 N + 308.599361) u v 3 + 495.005747 v 4 



(7) 



Pv 



v U- 



1 



N + 8 



'12u + 9v) + 



+ 92.444444 u v + 34.222222 



(3.407407 N + 54.814815) u 2 
1.251107 N 2 



(A^ + 8) 3 

41.853902 N - 469.333970) u 3 + 9 (0.248784 N + 136.511768) w 2 v 

1 



+ 957.781662 w v 2 + 255.929737 



(A^ + 



(0.574653 N 3 



- 0.267107 N 2 + 584.287672 N + 5032.692260) u 4 + 3 (0.057375 N 2 
+ 107.641680 N + 5989.283536) u 3 v + 3 (7321.464604 

- 16.494003 N) u 2 v 2 + 11856.956858 u v 3 + 2470.392521 v 4 



These equations are known to have four solutions corresponding to the 
trivial Gaussian, the Ising, the isotropic (Heisenberg) and the cubic fixed 
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points ||. The most intriguing of them are the isotropic and cubic ones. 
The one-loop approximation analysis of the eigenvalue exponents for these 
points yields the upper boundary value for the critical dimensionality N c , 
N C = A. 

To determine the lower boundary of N c one should attract the specific 
symmetry property of model (0), when N = 2 [ 26f . Namely, the transforma- 
tion of the field components 

<Pi ->• (<Pi + ¥2) , ¥2 -> -j= (<Pi - (P2) (9) 

combined with substitution of the quartic couplings 

3 , , 

u — > u + - v , v — > —v (10) 

does not change the structure of the initial Hamiltonian itself. As a result, 
the /3-functions (^), (§) should obey certain symmetry relations [p7| : 



/ 3 \ 3 

Pu[u+-V,-Vj = P U {U,V) + - P V {U,V) , 

f3 v (u+-v,-v^ = -(3 v (u,v) , (11) 



but their form remain unchanged. However, for N = 2 transformations (§) 
and (p!0| ) result in the relocation of the coupling constants values so that 
the cubic and Ising fixed points are transformed into each another at the 
3D RG flow diagram. Since the exact RG equations always have the Ising 
fixed point, which inevitably is the saddle-knot one, these equations should 
also have the cubic fixed point, which will be unstable. In this situation, the 
isotropic fixed point, again always existing in the exact RG equations, should 
be the stable knot only. Therefore, we conclude that the lower boundary of 
N c is not less than two. The real value of N c can be obtained only on the 
basis of analysis of the RG flow diagram structure, provided the /9-functions 
of the model are calculated in sufficiently high-order RG approximations and 
then processed by appropriate resummation techniques. 



3 Resummation, fixed points and stability 

It is well known that the field-theoretical RG series are divergent. The char- 
acter of their large-order behaviour is well established only for simple O(N)- 
symmetric models |28|, 0. The coefficients of the series at large k were 
shown to behave as c(—a) k k\k b , where the asymptotic parameters a, b and c 
are assumed to be calculated for each RG function. Knowledge of the exact 
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values of the asymptotic parameters in combination with the most powerful 
resummation procedure of the Borel transformation with a conformal map- 



ping, first proposed in Ref. [31], made it possible to develop the accomplished 
quantitative theory of critical behaviour of simple systems |3"5| . 

At the same time, the asymptotic nature of RG functions of anisotropic 
models is still unknown. Calculating the large-order asymptotic behaviour 
for the series in such models is a very difficult problem. That is why, lacking 
any information about the large-order behaviour either the simple Pade-Borel 
or Chisholm-Borel resummation procedures are used. The latter technique, 
however, possesses at least two inherent drawbacks. First, some ambiguity in 
the calculation of coefficients of denominators of the Chisholm approximants 



is unavoidable |34j . Second, the Chisholm-Borel procedure does not hold the 
specific symmetry properties of a model. At the same time, exploiting the 
Borel transformation in combination with the Pade or Chisholm approxi- 
mants shows that the results of calculation are very sensitive to the choice of 
the type of approximants. This may lead to estimates which do not provide 
reliable predictions even in the higher-loop RG approximations [ BfJ . Besides, 



in the framework of both schemes it is very difficult to determine any error 
bounds for the evaluated quantities. 

In this paper we attempt to overcome the outlined difficulties by apply- 
ing the PBL resummation method, generalized for the two coupling constant 
case, to processing the RG expansions (0) and @. This method, first intro- 



duced by Baker et al in |36| , turned out to be highly efficient when used to 
study the critical behaviour of the simple 0(A r )-symmetric models in 3D. The 
critical exponent estimates obtained within the framework of this technique 
are regarded nowadays as the most accurate values, as those of |37]]. 
We motivate our choice of the PBL resummation method with the following 
reasons. 

• 3D RG expansions for the /3-functions of the cubic model alternate in 
signs. Therefore, using PBL resummation technique is quite natural. 

• It can be expected that for complex models with more than one cou- 
pling constant, the asymptotics of the RG series at large orders will 
include a factor k\k b . The PBL resummation method removes diver- 
gences of this type. 

• The PBL resummation method allows one to determine the error bounds 
for the physical quantities to be calculated, in a natural way. 

The generalized PBL resummation procedure consists of the following. Let 
a physical quantity F(u, v) be represented by a double series 

F(u,v) = Y,fvU l v 3 , (12) 
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where coefficients /y ~ {i + + j) b at large orders — > oo), the ad- 
ditional parameter b being an arbitrary non-negative number to be defined 
below. Associated with the initial series (|12|) is the function 

oo 

F(u,v;b) = J e-H b B(ut,vt)dt (13) 
o 

The Borel-Leroy transform B(x, y) is the analytical continuation of its Taylor 
series 

which is absolutely convergent in a circle of non-zero radius. In order to 
calculate the integral in (|T3|) one should continue analytically B(x, y) for 
< x < oo and < y < oo. To this end, the rational Pade approximants 
[L/M] (x,y) are used. The Pade approximant method is determined in a 
conventional way |T3. Let us consider a "resolvent" series 

oo k e rfX„.k—l oo 

B(x, y, A) = V \ k V Sn , n = E , (15) 

where coefficients are uniform polynomials of fcth order in u and v. The 
sum of the series is then approximated by 

B(x,y) = [L/M]\ x=i . (16) 
The Pade approximants [L/M] in A are given by 

[L/M] ^ %§> • (17) 

where -Pl(A) and Qm(A) are polynomials of degrees L and M, respectively, 
with coefficients depending on x and y which should be determined from the 
conditions 

Qu {X)B{x, y- A) - P L (X) = 0(A L+M+1 ) , 

Qm (0) = 1. (18) 

Replacing variables x = ut and y = vt in the Pade approximants and then 
evaluating the Borel-Leroy integral 

oo 

F(u, v;b)= J e~H b [L/M} | x _dt (19) 
o 
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we obtain the approximate expressions for RG functions. 

Among Pade approximants the diagonal (L = M) or near-diagonal ones 



were proved to exhibit the best approximating properties [13]. However, as 
the degree of the denominator M increases, the number of possible poles of 
the approximant increases too. If some of the poles belong to the positive 
real semiaxis, the corresponding approximant should be rejected. Due to this 
the choice of "working" approximants, which might be used for analytical 
continuation of the Borel-Leroy image onto the complex cut plane, is largely 
limited. On the other hand, varying the free parameter b in the Borel- 
Leroy transformation (|l|) allows one to optimize the resummation procedure 
under the condition that the fastest convergence of the iteration process is 
achieved. So, taking into account the above-mentioned remarks, in order to 
find the locations of the fixed points we adopt the following scheme. For the 
fixed N, the /5-functions are resummed by virtue of transformation (|T3|) in 
the highest-loop orders by shifting the transformation parameter b. For an 
analytical continuation of the Borel-Leroy transforms B u (u, v), B v (u, v) over 
the cut plain the most appropriate Pade approximants [2/1], [3/1] and [2/2] 
are used. The locations of the fixed points are then determined for each b 
from the solution of the set of equations: Pl es (u c ,v c ) = 0, /3£ es (u c , v c ) = 0. 
The "true" locations are obtained by averaging over the values given by 
the approximants under the optimal value of the parameter b, at which the 
quantity [1 — jF L (u, v; b) / J-'i-iiu, v ; b)\ reaches its local minima. The quantity 
Tl{u, v; b) is evaluated for the L-partial sum of the series in equation ([19|), 
where L denotes the step of truncation of the series. 

The results of the computation of the cubic fixed point locations depend- 
ing on the parameter b are presented for the physically important case iV = 3 
in figure 1. Three curves correspond to the three Pade approximants. The 
parameter b shifts from to 3. As can be seen from the figure the opti- 
mal value of b is zero. At this point the numerical values of the cubic fixed 
point locations given by different approximants are the closest to each other. 
The result of computing the cubic fixed point locations for N = 3 are also 
presented in table 1. In the first three columns of the table the fixed point 
locations values found by means of the Pade approximants [2/1], [3/1], and 
[2/2] at b = are given. Averaging the results of processing over all of the ap- 
proximants under the optimal value of b gives the estimates presented in the 
fourth column of the table. We adopt these numbers as the final estimates 
for the cubic fixed point locations found within the four-loop approximation. 
As the degree of accuracy for these approximate values we take the maximum 
deviations of the average values of the fixed point locations from those given 
by the approximants at b = 0. 

One can observe, looking at figure 1, that the values of the cubic fixed 
point locations given by the symmetric approximant [2/2] depend weakly on 
the shift parameter b. Averaging over all the values given by this approximant 
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within the interval [0,3] results in the cubic fixed point locations estimates 
presented in the fifth column of table 1. The coordinates of the cubic fixed 
point found earlier on the basis of the three- and four-loop approximations, 
using the Chisholm-Borel resummation method, are presented in the table, 
for comparison. These numbers include the normalizing multiplier -g- needed 
to compare our /3-functions with those obtained in || jlOR . 

To verify the correctness of the chosen approach let us apply the above 
considered scheme to estimate the fixed point locations of the 0(3)-symmetric 
model for which the numerical results are well known. The six-loop 3D RG 
expansion for the /^-function of this model was reported in []3l], |36|. The 



PBL resummation of that series using eight types of Pade approximants [2/1], 
[3/1], [2/2], [4/1], [3/2], [5/1], [4/2] and [3/3] yields, after solving the equation 
f3 res (.9c) = 0, the picture displayed in figure 2. It is seen that the values of the 
isotropic fixed point location calculated in the highest RG orders with the 
help of the approximants [3/3], [4/2] and [3/2] are very weakly dependent 
on the parameter b varied within the interval < b < 15. The curves 
corresponding to these approximants are intersected at the point b = 4.5. 
Therefore b = 4.5 is the optimal value of the transformation parameter for 
which the fastest convergence of the iteration procedure is ensured. For 
b = 4.5 the central value estimate of the isotropic fixed point is g c = 1.392. 
The maximum deviation of the central value from the values given by some 
of the approximants [3/3], [4/2] and [3/2] at the point b = 10 is adopted 
approximately as an apparent accuracy of the calculation, A = 0.0013. Such 
a small error can be explained by the small dispersion of the curves within 
the range 5 < b < 10. So, the estimate g c = 1.3920 ± 0.0013 is in excellent 
agreement with those found more then 20 years ago in [31], 0] as well as with 



recent results of 32 



Within the framework of the four-loop approximation there are only three 
appropriate Pade approximants. Averaging the results of computing the 
isotropic fixed point location given by the approximants [2/1], [3/1] and 
[2/2] under the optimal value of the transformation parameter results in the 
estimate g c = 1.3925 ± 0.0070. The error was determined again through 
the maximum deviation of the central value from those given by each of 
the approximants at b = 0. It is seen that the four-loop estimate of the 
coordinate of the isotropic fixed point is in a good accordance with the best 
ones followed from the six-loop consideration. 

Note, however, that the coordinate of the 0(3)-symmetric fixed point 
calculated within the five- loop approximation does not approach the "exact" 
value. Namely, the PBL resummation procedure leads to the estimate g c = 
1.3947±0.0040. Although the error of the calculation became visibly smaller, 
the central value of the fixed point location stepped aside from the four- and 
six-loop ones. 

Thus, the fulfilled numerical analysis shows that the isotropic fixed point 
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location estimate obtained in the four-loop level occurs close to the six-loop 
value. One can expect, therefore, that in the case of the cubic model the 
fixed point locations u c = 1.3428 ± 0.0200, v c = 0.0815 ± 0.0300 (fourth 
column of table 1) will not be strongly distinguish from the "exact", say, the 
six-loop, values. The coordinates of the cubic fixed point for some N are 
presented in table 2. Our calculations show that for N = 3 the coordinates 
of the cubic fixed point practically do not differ from those of the Heisenberg 
one. However, with increasing N the cubic fixed point runs away from the 
isotropic point moving towards the Ising one. In the large N limit these 
two fixed points become close to each another so much that the influence of 
the 0(A")-symmetric invariant on the critical thermodynamics of the cubic 
model vanishes. This can be easily seen by applying the consideration to 
the one-loop solutions of the RG equations of model (p. Indeed, rescaling the 
coupling constants u ^ u/N, v ^ v/N in the initial Hamiltonian and taking 
then the limit iV — > oo one can see that the cubic fixed point approaches 
the Ising one asymptotically. So, the cubic model turns out to be split into 
N non-interacting Ising models, the critical behaviour of each of them will 
be determined by a set of the critical exponents renormalized according to 



Fisher 21 



Another way to determine the fixed point locations in fixed D is to con- 
struct the RG flows diagram of the model. If, at the flows diagram, there 
exists a fixed point of stable knot type, the trajectories originated from some 
point within the range of stability of the initial Hamiltonian would flow to- 
wards the knot. The region at the flow diagram where the trajectories are 
intersected provides the coordinates of the stable fixed point. Investigating 
the 3D RG flow diagram of model (|]) in the four-loop approximation we 
arrive at the conclusion that the cubic rather than the isotropic fixed point 
is absolutely stable for all N > 3. 

At the same time, the reliable prediction about the stability of the cubic 
fixed point for N > 3 can be made on the basis of calculating the eigenvalue 
exponents A's of the stability matrix 

Mij= U $L 

\ du dv / 

taken at u = u c and v = v c . If the real parts of both eigenvalues are negative, 
the fixed point is the stable knot in the (u,v) plane. If Ai, A2 have opposite 
signs, the point is of the "saddle-knot" type. 

To calculate the stability matrix eigenvalues of the cubic and isotropic 
fixed points we have chosen the following strategy. First, the derivatives 
of the /3-functions (0) and (|8|) are calculated, and the new RG expansions 
resummed by means of the PBL technique are substituted into the matrix 
Mij. The eigenvalue exponents of the matrix of derivatives My obtained in 
such a way are then evaluated under the optimal value of the transformation 
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parameter b. In figure 3 we present our numerical results for the series — ^ff 
and — for the physically interesting case N = 3. The curves correspond 
to the three types of the Pade approximants used within the four-loop ap- 
proximation. The crossing of the curves gives the optimal value of b at which 
we find ^\opt = -0.7536 and ^\ opt = -0.0331. Because the series -fjf 
and — ^r- turn out to be shorter by one order in comparison with — and 

au J L ou 

— ^ff, their resumming performed with the help of the approximant [2/1] 
only yields the monotonic dependence of the result of processing on the pa- 
rameter b. In this unfavorable situation, we take into account an additional 
Pade approximant [1/1] to optimize the iteration procedure. The results are 
plotted in figure 4. For the optimal values of b we obtain ^j^\ op t = —0.4566 
and Q^\opt = —0.0409. Straightforward calculation of the eigenvalues of the 
stability matrix Mij gives for the cubic fixed point the numbers provided 
in table 3. The eigenvalues of the isotropic fixed point as well as the anal- 



ogous numerical estimates obtained recently in |15| on the basis of using 
the five-loop e-expansions are presented for comparison therein. These esti- 
mates show that the cubic fixed point is absolutely stable in 3D for N = 3 
while the isotropic fixed point appears to be stable on the li-axis only. Our 
numerical results agree well with those obtained in flT5||. Unfortunately, at 



present we cannot indicate realistic error bounds in our calculation of the 
stability matrix eigenvalues. Nevertheless, a crude estimate can be done. In 
fact, if the model (|1]) is almost identical to some marginal system for which 
N c = 3 in 3D, the stability index A2 both for the isotropic and for the cu- 
bic fixed points should be equal to zero at iV = 3 and, consequently, the 
points should coincide. Therefore, as can be seen from the numbers given in 
table 3, the four- loop approximation predicts eigenvalues with error about 
0.01. Of course, dealing with the theory without a small parameter and the 
short perturbative expansions one would refer to such a level of accuracy 
as satisfactory. However numerically small errors may lead, sometimes, to 
qualitatively incorrect results |35| . 

Note, that although the recent high precision MC simulation using fi- 
nite size scaling techniques [[L6| predicts the stability of the isotropic rather 
than the cubic fixed point, the absolute value of the stability eigenvalue | A2 1 
obtained for the isotropic point turned out to be very small. This is in 
accordance with our estimate. 

Let us now calculate the critical dimensionality iV c of the order parameter 
field. The critical dimensionality is defined as a value of N at which the cubic 
fixed point coincides with the isotropic one. Equivalently, for N = N c the 
second eigenvalue of the stability matrix vanishes, A2 = 0. Studying 
carefully the 3D RG flow diagram of model (0) depending on the order of 
approximation for different Pade approximants we arrive at the conclusion 
that N c = 2.910 ± 0.035 and N c = 2.890 ± 0.020 within the three- and four- 
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loop approximations, respectively. The accuracy of the calculation of N c 
was determined through the evaluation of the stability matrix eigenvalues 
for different N from the interval of errors mentioned above. That value of 
iV = N c , above or below its central number, at which the second eigenvalue 
A2 was becoming non-zero, was taking for the upper or lower boundary of 
N c , respectively. 

It is worth comparing the four-loop estimate of N c just found with that 
which can be obtained within the e-expansion method. The five-loop e- 



expansion for N c has been calculated in | Tl| . The series turned out to be 
alternating in signs that allows one to resum it by means of the PBL tech- 
nique. To this end, we use again the most appropriate Pade approximants 
[2/1], [3/1] and [2/2] for analytical continuation of the Borel-Leroy transform 
for all < et < 00. Dependence of the results of processing of the critical 
dimensionality N c on the transformation parameter b is depicted in figure 5. 
The curves corresponding to the approximants are crossed at the point b ~ 1. 
The appropriate value of the critical dimensionality is N c = 2.894 ± 0.040. 
As an error of the calculation it is natural to assume the maximum scattering 
of numerical values given by the approximants at b = from that obtained 
at the crossing point of the curves. This estimate of N c is in excellent ac- 
cordance with the above found within the 3D RG approach. It also agrees 
well with the estimates obtained earlier on the basis of the different resum- 
mation technique || [UJ. So, both schemes, the RG technique directly in 
3D and the ^-expansion method, result in the same estimate of the critical 
dimensionality iV c = 2.89, thus implying that the cubic fixed point is stable 
in three dimensions for N > 3. This means that the critical behaviour of the 
magnetic phase transitions in crystals with cubic anisotropy should belong 
to the cubic rather than the isotropic universality class with a certain set 
of critical exponents. However, due to the obvious marginality of the model 
(N c ~ 3) and closeness of both (isotropic and cubic) fixed points on the 3D 
RG flow diagram for N = 3, the critical exponents values of the cubic point 
will be practically the same of the isotropic one. That is why the calcula- 
tion of the critical exponents in cubic magnets with N = 3 seems to be of 
academic interest only. 



Conclusion 

To summarize, the complete analysis of the global structure of RG flows of 
a model with two quartic coupling constants associated with isotropic and 
cubic interactions describing magnetic and structural phase transitions in 
a good number of real substances has been carried out within the massive 
field theory directly in three dimensions. Perturbative expansions for the (3- 
functions were deduced for generic N up to four-loop order. The fixed points 
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locations were found for iV > 3 by applying the generalized Pade-Borel-Leroy 
resummation technique. On the basis of comparative numerical analysis with 
the 0(iV)-symmetric models, the fixed point locations for which have been 
solidly established |32L B3J, we have made the assumption that the four-loop 
estimates of the cubic fixed point locations should not differ strongly from 
the "exact" values within the error bounds. 

The analysis of the eigenvalue exponents of the isotropic and cubic fixed 
points fulfilled for the physically significant case N = 3 has shown that 
the cubic rather than the isotropic fixed point is absolutely stable in 3D. 
The eigenvalues estimates (see table 3) were found to agree well with those 
calculated on the basis of exploiting the five-loop e-expansions combined with 
a careful resummation procedure [15| . Our results agree numerically with the 
recent high precision MC estimates [16|, in spite of the latter predicting the 



stability of the isotropic rather than the cubic fixed point. 

The critical dimensionality iV c of the order parameter, at which the topol- 
ogy of the flow diagram changes, has been estimated by the two different 
methods: (a) by resumming the four-loop RG expansions for the /3-functions 
in 3D and (b) by resumming the five- loop e-expansion for N c at e — 1. The 
numerical estimates N c = 2.89 ±0.02 and N c = 2.894 ±0.040 obtained are in 
a good agreement with the earlier results |Tl]] and confirm the conclusion 
about the stability of the cubic fixed point for N > 3. Consequently, the 
magnetic and structural phase transitions in three-dimensional anisotropic 
crystals with cubic symmetry are of second order and their critical thermo- 
dynamics should be governed by the cubic fixed point with a certain set 
of critical exponents. Unfortunately, the cubic universality class is not eas- 
ily distinguished experimentally from the isotropic one, due to the obvious 
marginality of the problem, iV c ~ 3. 
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Figure Captions 



Fig. 1. Curves demonstrating the dependence of the results of calculating 
the cubic fixed point locations on the transformation parameter b for N = 3. 
The upper curve (O) corresponds to the [2/1] approximant, while the middle 
(A) and lower (□) curves correspond to the [2/2] and [3/1] approximants, 
respectively. 

Fig. 2. The results of computation of the 0(3)-symmetric fixed point loca- 
tions from the three- to the six-loop approximations obtained on the basis of 
the PBL resummation method with eight types of the approximants: [2/1] - 
O, [3/1] - □, [2/2] - A, [4/1] - full O, [3/2] - full A, [5/1] - full o, [4/2] - o, 
[3/3]- x. 

Fig. 3. Graphs of dependence of the results of processing of the series a) 

— , b) —^§f on the parameter b, N = 3. The curves are given in the same 
notations as in the previous figures. 

Fig. 4. Graphs of dependence of the results of processing of the series a) 

— ^o^i b) — ^ on the parameter 6, N — 3. For the curves corresponding to 
the approximants [1/1] and [2/1] the notation O and □ are used, respectively. 

Fig. 5. Dependence of the results of processing of the e-series for the critical 
dimensionality N c on the transformation parameter b. 
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Table 1: Coordinates of the cubic fixed point of RG equations for N = 3 
found under the optimal value of the transformation parameter b = 0. 





[2/1] 


[3/1] 


[2/2] 


Average value 


Average over [2/2] 




1.3536 


1.3338 


1.3410 


1.3428 ± 0.0200 
1.3480°, 1.3357 6 


1.3425 


V c 


0.0526 


0.1026 


0.0894 


0.0815 ± 0.0300 
0.0904°, 0.0906 6 


0.0937 



a Quoted from Ref. g 
b Quoted from Ref. |§ 



Table 2: Coordinates of the cubic fixed point of RG equations for some N 
found under the optimal value of the transformation parameter b within the 
four-loop approximation. The average values of the coordinates calculated 
over the most stable Pade approximant [2/2] are also presented for compar- 
ison. 



N 


4 


5 


6 


7 


8 


9 


10 


U c 


0.9055 


0.6980 


0.5807 


0.5060 


0.4544 


0.4168 


0.3881 


V c 


0.8167 


1.2361 


1.5386 


1.7874 


2.0076 


2.2108 


2.4032 


u c [2/2] 


0.8981 


0.6886 


0.5708 


0.4962 


0.4448 


0.4074 


0.3789 


v c [2/2] 


0.8380 


1.2608 


1.5649 


1.8148 


2.0359 


2.2400 


2.4333 



Table 3: Four-loop eigenvalue exponents estimates for the cubic (CFP) and 
isotropic (IFP) fixed points found for N = 3 under the optimal value of the 
transformation parameter b. 





CFP 


CFP, Ref. [15] 


IFP 


IFP, Ref. P| 


Ai 


-0.7786 


-0.7648 


-0.7791 


-0.7640 


A 2 


-0.0081 


-0.0085 


0.0077 


0.0089 
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